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ABSTRACT 


The Aero-As trodynamics Research Review No. 2 (NASA TM X-53295) con- 
tains a succinct account of Newtonian aerodynamics for general analytical 
surfaces. Much of the detail passed over there when dealing with appli- 
cations, together with a more thorough exposition of the general approach, 
is given in the present report which, it is hoped, will be useful to those 
who wish to actually work with Newtonian's method. It ought to serve well 
in the high-Mach number, low-density phase of re-entry flight. The 
lemniscatic body chosen in the Research Review paper as an example for 
calculating blunt-nose forces is replaced here by a circular cone with 
spherical base cap, since this configuration seems to command more immedi- 
ate interest. 
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LIST OF SYMBOLS (Partial) 


A, General : Vector quantities are denoted by underlining, 

except the unit vectors 


i, j, k 

on the axes of a Cartesian rectangular system. 

CX: overall angle of attack 

C^ oc : angle of attack at an elemental surface 


a % : its complement 

n: unit vector in direction of the surface interior 

normal (components, n x , n 2 , n 3 ) 

v: unit vector in flow direction (components a x > Q^, C%) 

Poo> Vqq 5 density, velocity, dynamic pressure in undisturbed 

flow 

Cp: local pressure coefficient 

a, t: coordinates in a surface system 

dS: elemental surface 

d 2 JP: elemental force acting on dS 

P: total force (components X, Y, Z) 

0: origin and point of reference for force moments 

xz radius vector from 0 to dS 

M: moment integral of elemental forces 

M*: moment of total force 

r*: its arm 

s t : tangential shadow curve 

s c : cast shadow curve 

iv 



B. Elliptical Cone ; Origin at tip; flow toward it. 


a: length unit on x-axis (which is taken as cone axis) 

a: surface coordinate, here also serving as multiplica- 

tor to express abscissas x in terms of a 

b, c: elliptical semi-axes at a = 1; b > c 

B, C: semi-axes of cone base (where a = o^) 

e: numerical excentricity of the elliptical cross 

sections 

co: angle made by top meridian and cone axis (tip 

angle) 

q = e 2 cos 2 w 

JU 

x : abscissa of centroid if placed on cone axis. 


C. Coaxial Cylindrical Afterbody : Origin at cone tip (as above). 

S', semi-axes of its cross section, parallel to B and C; 

B, C ^ B, C 

A = (ff/B) 2 

a, t: surface coordinates on cylinder 

x^ = aa^: cone length 


x c = a o c : total length of composite body 

h = a(a c - a^): cylinder length 

a = 'i)t(t): equation of cast shadow curve on cylinder 

t = 0: tangential shadow line on cylinder 

R: cross sectional radius if cylinder is circular; 

R ^ C 


0 = R/B 


v 



Z: force acting on cylinder 

3: its moment with respect to cone tip 

x*: abscissa of cylinder centroid 

x' : centroid abscissa of composite body (cone and 

° afterbody) 


D. 


Biparabolic Conoid : Origin at tip; flow toward it. 

a: parabola parameter, used as unit length 

go: angle made by top meridian and conoid axis (tip 

angle) 

/V A 


B, C: 

% = 2a °b : 
X: 


"semi-axes” of biparabolic base (Figure 4) 
conoid length 

wedge angle made by the x,y-plane and any plane 
through y-axis 


m = \I (x b / a) 


sin 2o) 


(1 + m 2 ) arc tan m - m 

A - m 3 ' 

K = n/i + (xb/Bp 
B: see Elliptical Cone 

x' c : abscissa of centroid if placed on conoid axis. 

E . Circular Cone with Spherical Sections Capping Its Base : Or ig in 

at sphere nose; flow toward it. 

B, R: cone base and sphere radii 

a; cone length 

xu: distance of cone base from origin 


vi 



P : 


surface coordinates on sphere 
surface coordinates on cone 


a, p: 

p*: defines tangential shadow meridian on cone 

X, Y: y-components of the total forces acting on spherical 

cap and cone, respectively 

x*: centroid abscissa of cone alone 

x tot : centroid abscissa of composite body (capped cone) 

The Modified Approach : a « 0°. 

0 ): local angle of attack at stagnation point 

it 

C p : pressure coefficient at stagnation point 

C p : pressure coefficient elsewhere 

M : Mach number in undisturbed flow 

CO 

y: ratio of specific heats in the gas 

q q : angle between attached shock and axis of a pointed 

body of revolution (shock angle) 


K = M a> 


K = M a 
s 00 s 


vii 
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NEWTONIAN AERODYNAMICS FOR GENERAL BODY SHAPES 
WITH SEVERAL APPLICATIONS 

SUMMARY 


In the hypersonic regime the Newtonian flow model, especially in 
its modified form, has been known for some time to produce satisfactory 
results regarding a number of basic body shapes. It may reasonably be 
expected to also work well with more complicated body geometries as 
presented, e.g., by re-entry vehicles. Because of its simplicity, the 
components of the aerodynamic force and the location of the centroid can 
be calculated without the elaborate effort otherwise required in gas 
dynamics. The mathematics pertaining to the unmodified approach have 
been developed here for a general surface given in analytic terms. Appli- 
cations have been worked out to the elliptic cone, first without, then 
with an attached rear cylinder, to a conoid of biparabolic cross sections, 
and to a blunt body resembling the Apollo capsule. The modification merely 
amounts to improving a constant; it has been described in the last section 
where a survey of results as compared with known data is also given. 


I. INTRODUCTION 


The Newtonian concept substitutes for the flowing liquid an assemblage 
(point mass array) of many minute inelastic particles all moving with equal 
and parallel velocities. On contact with a material object such particles 
transfer their momenta components normal to the surface while they retain 
their tangential momenta carrying them off without further effect unless 
there is a secondary impingement by the deflected stream. There is no 
action on the body except by direct hits; no pressure exists on surface 
elements shielded by upstream elements (contrary to later experimental 
evidence). The formation of a shock, which of course was unknown in 
Newton's time, is even now not contemplated in simple impact theory. 

This theoretical picture of a flowing liquid and of its actions on 
body surfaces had soon to be abandoned in hydrodynamics. It is also 
inadequate with gas flows, both subsonic and supersonic. But in the 
hypersonic regime, the expressions derived from it for the surface pres- 
sure distribution and concomitant aerodynamic data yield values often 
surprisingly close to those observed experimentally or obtained by more 
exact theoretical means. In recent times Newton's method is freely used 


in dealing with such flows; it offers closed-form results available in 
a relatively easy manner. It has been found wanting, however, with sur- 
faces concave to flow direction and in gases with exceptionally low ratio 
of specific heats approaching unity (for explanation and examples, see 
Reference 1, p. 125-128). Applications seem best restricted therefore to 
bodies turning convex surface parts toward an onrushing gas of ordinary 
description. The simple theory can often be improved upon by introducing 
a correction to the local pressure coefficient based on shock transition 
relations (modified method). A dependence on Mach number, absent so far, 
then appears. 

A hypersonic regime should prevail for some time when a space 
vehicle, approaching the earth with quasi-cosmic velocity, moves in the 
uppermost regions of the continuum atmosphere. One would like to pre- 
dict theoretically the lift and drag the body will experience in these 
regions, together with the stability behavior as depending on the loca- 
tion of the center of pressure. 

However, the body shapes considered in the literature are mostly of 
the simple kind, such as circular cones and cylinders, spheres, diamond- 
shaped wing profiles, blunted cones and capped cylinders. There appears 
to be a need to systematize the method for application to more complex 
structures. In the following, the Newtonian expressions will be derived 
for an unspecified general surface which is allowed to be composite. 

The stipulation is made that all its parts can be described by analytic 
equations. Surface zones hit by secondary incidence are not permitted 
to exist. Applications will be made to the elliptical cone (with and 
without cylindrical appendage), to a conoid of biparabolic cross sections, 
and to a blunt-nosed body of the Apollo capsule types. 


II. IMPACT FLOW MATHEMATICS IN GENERAL 


The angle of attack, a, is commonly understood as the angle made by 
the direction of the uniform flow with a line chosen within the body, 
usually the body’s axis if such an axis can be defined. In conventional 
aerodynamics there is rarely occasion to consider other angles relative 
to flow direction, since it affects the pressure distribution only 
indirectly by way of boundary conditions. But in impact flow where every 
elemental surface, if hit at all, is hit from a well-defined direction 
bearing on the elemented force transferred, a local angle of attack must 
be defined which will depend on the position of the local surface element 
as indicated by its interior normal. Let the latter’s direction be given 
by the unit vector 


n = n x ± + n 2 j + n 3 k 


( 1 ) 


2 



where i, j and k are unit vectors on the axes of the Cartesian system in 
which the body's surface is described. The unit vector, _v, in flow direc- 
tion will usually be linked to the position of the axis; its expression 
will then contain certain trigonometric functions of a which cannot be 
obtained before a particular body-flow configuration is considered. At 
this station we therefore write, in general, 


v_ « a x i + a 2 j 4- a 3 k. 

The two vectors n and v_ define the local plane 
normal to the elemental surface and intersects 
in a straight line, t. The angle t makes with 
general, considered the local angle of attack, 
handle, however, is its complement 

= f - “loc < 3 > 


( 2 ) 

of incidence which is 
with its tangential plane 
the direction v is, in 
°loc* More convenient to 


since as the acute angle between flow direction and the interior 

normal ", is easily obtainable from the scalar product 


cos a' 


n • v 


I 


n.a. . 


1 1 


(4) 


In this way one eliminates the somewhat bothersome task of determining 
the line t which, except for C^ oc> has little interest, whereas n gives 
the direction of the local force. To be sure, the tangential momentum 
of the striking particle is carried away in the direction of t, but that 
momentum component is ineffective, secondary hits not being admitted. 

Aside from terms small of second order, the local pressure coef- 
ficient will be that of a flat plate under the angle of attack “Loc- 
Newton's second law, if applied to normal momentum transfer, gives the 
express ion 


Cp 2 s in C^i oc 


(5) 


"The interior normal, if defined as the axis of the elemental surface, 
would give of the meaning of local angle of attack in agreement with the 
definition of the overall angles. 
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which, through relations (3) and (4) goes into 


Cp = 2(v * n) 2 . 


( 6 ) 


Expression (5) is well known (Reference 2, p. 6); it uses the continuum 
notion of density to describe the mass of the particles per unit volume 
of the incoming flow. The elemental force, as usual, is obtained from 
the coefficient Cp by multiplying it with the local surface element, dS, 
and with the dynamic pressure. 


q 


00 


V 2 


2 


(7) 


of the undisturbed flow which is taken as known. Thus, 


d 2 P = 2q (v • n) 2 dS n. 

— x oo — 1 — — 


( 8 ) 


The local force is small of second order and points into the direction n. 
The determination of n and dS calls for the methods developed in the dif- 
ferential geometry of surfaces. Here, lucid and symmetric formulations 
emerge when the surface is given a point-wise representation 



(9) 


The variables a and t will move within certain "natural” intervals in 
order for the Cartesian triplets x, y, z to exactly embrace the surface 
points. Often, a and t will be found by geometric reasoning which may 
lead to a clearer visualization of surface features and, by the way, may 
suggest several attractive ways of introducing these variables. On the 
other hand, one might be quite formal about it, putting x = a, y = x, 
z = f(a,T). 

Representation (9) offers the additional advantage that, on assigning 
parametric values to t or to a, it defines two sets of parametric curves 
on the surface 

x x = x(ct,t*), y x = y(a,T*), 
x 2 = x(ct*,t), y 2 = y(c*,x). 


Zi = z (a, t*) 
z 2 = z (a’% t) 


4 


( 10 ) 



which constitute a curvilinear (not necessarily orthogonal) system of 
coordinates on it. The lines (10) are individualized by the choice of 
T* and a* at which point they intersect. In general, they must be 
expected to be of double curvature (spatial curves). However, they can 
often be taken as meridians and cross-sectional peripheries, in which 
case both sets are planar. If in doubt one may consult the determinant 


X* 

y' 

z' 



x" 

y" 

z" 

= A 

(11) 

x". 

ylll 

z"' 




where the primes denote differentiation with respect to the one quantity, 
a or t, which is variable in the set considered. If A = 0, that set con- 
sists of plane curves. It is seen that if, e.g., the surface representa- 
tion (9) depends in the second degree at most on a, the curves t* = const, 
will be planar. One can then use one of the planes for a two-dimensional 
representation of the curve it contains and thus make the latter 1 s shape 
and course on the surface more clearly understood. 


The surface functions (9) depend on both the variables o and t. 
Partial derivatives with respect to them will be indicated by correspond- 
ing subscripts. 

Regarding first the vector n, its components follow from the set (9) 
as 


1 

z a 

i 

z a x cr 

1 

z a y CT 

ni - ±N 

y t z t 

5 112 " ±N 

Z T X T 

’ n3 _ ±N 

Z T 


where \ (12) 




2 


2 


2 

/ 

Vo z a 

+ 

z o x a 

+ 

X a ya 


J 

y T 

z T x T 

yx 
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The sign appropriate to the interior normal can be fixed without difficulty 
once the functions (9) are known for a particular surface. 

Secondly, the elemental surface is given by 

dS = N da dr (13) 

where da, dx should be introduced as positive increments (the later inte- 
grations must be performed over intervals in which a and t increase). 

The second order force differential (9) may now be written as 


d 2 P = 2q oo (n 1 Q! 1 + n^X 2 + n 3 a 3 ) 2 N(in x + jn 2 + kn 3 ) 


= i d 2 X + j d 2 Y + k d 2 Z 


(14) 


where the components X, Y, Z of the total force P are often called the 
"chordwise," "lateral," "normal" forces, respectively. These terms 
imply that the chord is taken as the x-axis and that v_ has no component 
in y-direction (rather than z-direction) . They must be exchanged with 
different correlations, and may even lose meaning altogether with highly 
irregular bodies. But in any case the determination of P is reduced to 
a double quadrature which can always be carried out either analytically, 
or, failing that, numerically. 

In addition to the resulting force P there will be a resulting 
moment of the elemental forces: 


M = 


[r X d 2 P] 


(15) 


where 


r = ix(cr, t) + jy(a,x) + kz(cr,x) 

is the lever arm from the origin (serving as the reference point) to the 
elemental surface, dS(o,x) where d 2 P is attacking. 
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From the condition that M should be equal to the moment of the 
resulting force. 


M* = [r* X P], 


( 16 ) 


the arm 


r* = ix* + jy* + kz* 


of the resulting force can be calculated. However, there are infinitely 
many of such arms, since a force may be moved along its line of attack 
without changing either the translatory or the rotatory effect. Mathe- 
matically, one of the r^- components remains indeterminate. An additional 
condition 


f(x*, y*. z*) - 0 (17) 

may be arbitrarily set up which in symmetrical configurations usually 
follows from the desire to have the point of r* (the centroid) on the 
body axis. 

Integration Limits 

While the foregoing formulas are perfectly general and directly 
applicable to any analytic surface (or surface part), the determination 
and proper treatment of shielded areas often require painstaking detail 
work, especially if there are several such areas (which may or’ may not 
overlap). The problem remains simple when the "shadow" is a single 
point or an open line (enclosed area zero), because then the reduction 
in force is a zero quantity, and the natural a- and t- intervals can still 
be used when integrating. 

Two types of shielded zones can be distinguished: those created by 

the bulk of the body as it opposes itself to the stream (tangential shadow 
curve), and those caused by the shadow cylinders through this curve and 
through sharp rims or edges that may exist on the surface (cast shadow). 

With the first type, the grazing flow vector will reside in a plane 
tangential to the surface and thus be perpendicular to the interior normal 
at the point of contact. It then has to satisfy the condition 

v • n = 0. (18) 
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This relationship of a and x, if written as x = <p(a) and introduced into 
the function (9), defines the points of the tangential shadow curve, s t , 
in terms of a. 

The natural limits of the a- and x- intervals, if traced in a (cr,x)- 
diagram delineate a usually rectangular image of the surface. The curve 
x = cp(cr) cuts off impact-free portions from this rectangle which must be 
excluded in integrating the force and moment differentials and thus changes 
the integration limits. Sometimes equation (18) will have the simple solu- 
tions a = const, and/or x = const., so that the integration area remains 
a narrowed-down rectangle. In the more complex cases involving a curved 
shadow boundary the integration will often have to be carried out numeri- 
cally. 

The root line of a shadow cylinder, as any spatial curve, can be 
described through the variation of a single variable, g: 

x c = y c = fsa) > Z c = f3a) - 

Care should be taken to determine the interval in which g may move in 
order to account for that length of the line that is exposed to the 
particle stream. Xf the curve is part of the surface (9) , as with the 
rim of truncated bodies, the variable g can be identified with either o 
or x (fin edges, e.g., would not usually permit this). 

The straight lines 


x - x c = y - y c = z " z c 

a x oc z a 3 


(19) 


in flow direction define the curved surface of the shadow cylinder; their 
intersections with other portions of the body surface, if they exist, 
form an additional shadow boundary, Sq, which, in handling the integra- 
tion limits, must be considered together with any shadow boundary, s^., 
consequent to tangential incidence and determined by equation (18). An 
example of this is found in Section IV. 
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III. THE ELLIPTICAL CONE IN SIMPLE IMPACT FLOW 


With the tip of the cone at the origin and its axis coinciding with 
the x-axis (Figure 1), its equation can be written as 


y- z 2 x 2 

-f — ^ £ 

b a 


( 20 ) 


implying that the cross section at any station x = const, is always an 
ellipse with axis ratio b/c. The lengths of the half axes at the base 
are denoted by B and C. Surface points with x < 0 will not be considered. 


Z 



Figure 1. 


The Elliptical Cone 
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One verifies easily that the coordinate triplets connected by equa- 
tions (20) can be represented as 


x = a cr I 

y = ba cos r > . (21) 

z = ca sin t J 

The variable a here is a pure number that counts the distance x from the 
origin in terms of the unit a. If a * identifies the distance of the 

base (xb = acr^ being the cone length), the variable o will range in 

0 ^ a ^ a b 

which is its natural interval. The parametric lines a = o* (second 
family is the set (10)) are cut out by planes parallel to the (y,z)- 
plane; they are the cross sectional ellipses of the cone. Condition (11) 
for planar curves is satisfied, since the differential quotients of x 
with respect to t are all zero. 

The surface lines r = x* (first family) are also planar (x, y, z 
depend linearly on a). They pass through the origin, since a can become 
zero. Their planes 


cy sin = bz cos t'% 


(22a) 


being parallel to the x-axis, therefore intersect on it. Any one of these 
meridional planes defines two generating lines, members of the first family. 
If it makes the angle, p, with the (x, y)-plane (Figure 1), it can also be 
described by the equation 


z = y tan p. 


(22b) 


It follows that 


b tan p = c tan t 


( 23 ) 


where t has been written for t*, since x“ can be any of the angles t. 
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To clarify the meaning of t geometrically, consider Figure 2, which 
uses the cross sectional area at the distance a = 1. The plane (22b) 
intersects with the elliptical circumference at point R # The point R 
is on a circle of radius b concentric with the ellipse and has the 
y-coordinate in common with R(y = y). To show that the angle t, as 
indicated in Figure 2, bears the relationship (23) with p we note that 
the condition y = y, if introduced into the equations of the circle and 
the ellipse, requires that 


z : b = z : c. 


But z = b sin t, so that z = c sin t. Furthermore, y = y = b cos t. 


Z 



Figure 2. Relationship of t and p 
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Finally, 


z c 

— = tan p = T" tan t 
Y b 

which is the relation (23). 

The meaning of the angle t now being seen, its range when R moves 
round the full elliptical circumference may be given as 


-It = T ^ K. 


(24a) 


If for simplicity we assume that the flow is parallel to the (z, x)-plane. 


v = i cos a + k sin a. 


(25) 


the incidence will be from below and behind as long as 0 ^ t S jt/2. The 
right half only of the cone needs to be considered in that case, provided 
that the chordwise and normal components are doubled (the lateral com- 
ponents cancel, Y = 0) . The t- interval is then narrowed down to 




(24b) 


The components of the interior normals are found from the general 
expressions (10) as applied to the surface (21): 


n n = 


bca 
N * 


no ~ ~ 


ac q cos t 
N 


n-a: — “ 


ab a si n T 
N 




N = + a \] b 2 c 2 + c 2 a 2 cos 2 t + a 2 b 2 sin 2 T 


(26) 


The signs have been chosen such that the x-component is always positive, 
as it evidently must be (Figure 1). 
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The shadow of the cone f s base rim does not strike the surface and 
may be disregarded. But a tangential shadow boundary may exist. Con- 
dition (18) can be put into the form 

sin t = — cotg a. (27) 


Since the right side is constant, the shadow line will coincide with the 
generatrix t = T t . Now by Figure 1, 


tan 



T 



Consequently, there is no real solution for if OC < w; no shadow 
develops, and the cone is impinged on in its entirety. The T- integra- 
tion will have to extend over the natural interval (24b). However, 
when OC > w, this region shrinks to 



T ^ 


T 


t 


with x t tc/ 2 if oc — > u). The part of the upper surface half lying between 
the generatrices t = T t and t = rt/ 2 is then shielded from the flow; it 
grows larger with increasing values of OC and comprises the entire top of 
the right half-cone if OC = 90° (x t = 0). 

The a- interval is evidently not affected by the presence of the 
shadow line Sf 

In writing down the expressions for the two non-zero force components 
it is convenient to use the abbreviation 


where 


= e c 


cos 


€ 



(b > c) 
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is the common numerical eccentricity of the ellipses. It then follows 
from the expressions (14), (25), (26) that 


, (sin w cos O' - cos w sin a sm t) 2 , 
d a = -4q bca da = — 2 — - — — dx 

00 1 - q COS T 


d 2 Z = -4q aba da 

^co 


(sin co cos a - cos cj sin PC sin x)‘ 


q cos^x 


sin t dx. 


On integrating over the intervals 0 ^ a ^ a^ and — at/ 2 ^ x ^ T t ,the force 
coefficients, if referred to the base area 


it be a 2 = jtBC , 


are found as 


2 sin 2 q J^r 
* e 2 (2 ' 't 


C„ = ~ — 72^ 1*$ + T. + \Fq tg OJ cotg a In 


1 + \Tq cos t 
1 - cos t 


+ ^ 1 1 [ tg^co cotg^a 

-fi - c 


- ±=-*)h + 


tan t 


arc tan 


^TT 


. (28) 


C = — — — jcos T + ^ (tg 2 w COtg^a 

Z it 6 tan to t 2 V 


- \ 1 + n/cT cos t 

In £ 

^ ' 1 - */q” cos t 


4- 2 tg w cotg OC 


~ + t - n/ 1 - q ( — + arc tan 


2 t 


tan t. 
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TC Ti 

The principal value (in <— ^ of the arc tan-function must be taken 

here. The numerical value of T t follows from relation (27). Expressions 
(28) need be used only with q > to. They simplify considerably when there 
is no tangential shadow zone on the surface (q ^ go; T t = it/ 2) and can then 
be set into the form 


c x = 2 


H 1 


- Q . p o COS 

— ^ s in q cos w H 


nT 


~q sin 2 cal 

1- q J 


> a ^ oo. 


2/v 1 - \ll - q 


C = 2 sin 2a cos^cj 
Z q 


(29) 


J 


The known formulas for the circular cone (Reference 2, p. 80-84) are 
obtained from expressions (28) and (29) by letting q approach zero. 


On account of the symmetrical configuration, the total force 
P = iX + kZ is in the (z,x) -plane and can be split there in two com- 
ponents, one in the direction of incidence (drag), the other perpendicular 
to it (lift), the coefficients becoming 


C = C cos a + C,, sin a 
D X Z 


C = - C sin q + C cos q 

L X Lt 


(30) 


There is no lift with q = 0. 

By somewhat laborious integrations it can be shown that the moment 
(15), needed for the determination of the centroid, has neither an i- 
nor a k-component. From the symmetry of the flow-body configuration, one 
would indeed expect that the moment should seek to turn the body about 
an y-axis parallel. It attains the form 


M = | q m j b(a 2 + c 2 ) a 2 F(x t ) 


15 


where the integral 



(sin a) cos a - cos 00 sin oc sin t) 2 
1 - q cos^t 


sin t d r 


also appears in the Z-component of the total force which may be written 
as 


Z = - 2q ab O'f F(t ). 

00 b t 


The moment (16) of the total force emerges as 


M" = 


J 

y* 


k 

jl 

z'' 

z 


= y w (iZ - kX) - j(x*Z - z*X) 


Requiring M = M one sees that y~ must be zero. This is understandable, 
since the force P is in the (z,x) -plane, its arm r* (relative to 0) must 
have the same property. Moreover, if we let the force attack at a point 
of the body- (x-)-axis, the arm will be part of the latter 


z 


= 0 . 


This is the form taken here by the condition (17). The equality of 
moments then leads to 


x 


lal _ 2 . c= . 2 % 

-Z 3 a °b ^ a 2 ^ 3 cos^co 


(31) 
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The location of the centroid depends only on the length, x^, of the cone 
and on the half-angle to. With these parameters preserved, it does not 
vary if the angle of attack or the elliptical shape, or both, are changed. 
The formula is known for circular cones, but often derived employing the 
normal forces alone* and then does not contain the denominator cos 2 oo, 
which indeed is negligible for small values of to. Yet it will be noted 
that, with a 35-degree cone, the centroid is already close to the base 
and moves out of the body if to is increased further. Since the mass 
3 

center is at ^ x^, slender cones are not stable. Stability begins at 
to ^ 20 degrees . 

The restriction to the symmetric case maintained so far is tantamount 
in real flight to permitting pitch angles only. If it is desired to con- 
sider yawed and rolled positions as well, one may use the same body-fixed 
(x,y,z)-system as before. The unit vector, v_, of incidence, however, will 
have the general form (2) where the direction cosines 0Cj_ are determined by 
the instantaneous flight direction. The simple shadow boundary condition 
(27) will be replaced by 


X -I a t ■ °- 

and the discussion of integration limits will become less straightforward; 
both lower and upper boundaries of t might be affected. The integrations 
can still be carried out analytically, as the force differential assumes 
the form 


d 2 P = 2abcq 


(Oh. + 

K a b 


cos X 




sin x) c 


_1 

a 5 


, cos 2 x . s in 2 x 


1 COS T 


k s in x 


dx a da. 


Closed expressions for the force components can be obtained from a table 
of indeterminate integrals. The Y-component is no longer identically 
zero, and the full range (24a) for x must be retained when seeking out 
the integration limits. 


/C See Ref. 2, p. 84. Mr. E. Lins ley seems to have been among the first to 
insist that the ch'ordwise forces can make a significant contribution 
to the moment M of a circular cone. 
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IV. CYLINDRICAL AFTERBODY ATTACHED TO ELLIPTICAL CONE 


The surface (21) becomes composite if, e.g., a coaxial elliptical 
cylinder is attached to the cone base. Let its cross-sectional semi- 
axes, £ and C, be parallel to, but neither larger than nor necessarily 
proportional to B and C (Section III, Fig. 1). With larger-base or 
rotated cylinders, the description complicates without calling for 
materially differing concepts. 

Since the cone and cylinder axes coincide, the same length variable, 
a, and the same angle p may be used for both surfaces. But, excepting 
the case where B : C = B : C, the variable T must be replaced by a vari- 
able t which will be connected with p by the equation 


B tan p = C tan t 


analogous to the elliptic relation (23). The points on the cylindrical 
surface 


p 2 

+ 2 - = 


1 


may be represented in the form 



(32) 


The same two sets of planes as were used with the cone determine the 
parametric lines t = T* and a = a* on the cylinder. 


The variable a ranges in the interval 


a, ^ a ^ a 

b c 

if aa c is the total length of the composite body. That of the cylinder 
alone will be 


h = a(a c - a b ) = * c - V 
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If again the symmetr ical flow vector (25) prevails, it suffices to 
study the right half of the cylinder. The natural interval in which the 
variable t moves then reduces to 


Both of the intervals will be curtailed by shadow boundaries (except with 
OL = n/2, when that of o remains intact). Before establishing the inte- 
gration limits, the force differential will be set up. 

From the general expressions (12), the interior normals are found as 


= 0 


n 2 = - 


a C cos 


N 


n 3 ~ 


a B sin t 
N 


where 


N = a 


n/b 2 


s in 2 T 


+ CT‘ 2 


cos 2 t 


On the right half of the cylinder, the component n 2 always points into 
the negative y-direction (n 2 < 0). 

It follows that 


v • n 


a B sin T sin a 
N 


(33) 


As with the cone, the lateral forces cancel (Y = 0). The chordwise 
forces are also zero since the elemental forces are all perpendicular to 
the x-axis; indeed, n 1 = 0. Expression (14) reduces to 


« 3 ** 

d 2 Z = - 4q a B 3 s in 2 a — do dr (34) 

B 2 sin 2 T + C 2 cos 2 r 

where the right side has been doubled to compensate for the use planned 
of only 1/2 the t- integration interval. The cylinder is affected by the 
normal forces only. 
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The tangential shadow line on the cylinder obviously is the genera- 
trix t = 0 (as can also be seen from the tangential condition • n = 0 
which here requires = 0) • The integration interval for x therefore is, 
at most 

- f s s °> 

it will be even narrower than that when the cast shadow reaches out beyond 
the cylinder length. 

Shielding through the cone f s tangential shadow cylinder does not 
occur as one infers from inspection. The somewhat lengthy analytic proof 
is omitted here. However, if B < B, or C* < C, the cone's base rim will 
cause an edge shadow curve to appear on the cylinder's lower surface, 
cutting off it an unimpinged zone. 

The coordinates of points on the elliptic rim line may be written as 


x r = *b 


y = B cos t 
r } 


z = C sin t 
r 


(35) 


The flow vectors passing through thejipper points of the base periphery 
miss the cylinder, and, except when B = B, even the shadow cast by an 
outer segment of the lower periphery will not strike it. For integration 
purposes the variable t thus will move in 


T ^ T ^ T 
2 max 


where the value of x max is zero or negative. It is associated with the 
particular generatrix of the shadow cylinder that is tangent to the 
material cylinder; i.e., that touches it somewhere along the line t = 0. 

The shadow generatrices (19) here have the form 

x - x z - z 

r _ r 

^ cos a sin a # 
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If the cylinder triplets (32) and the rim triplets (35) are introduced, 
these relations go into 


B cos t = B cos T 

(x - x^) sin O' = (C sin T - C sin t) cos OL 

where the pairs x^,, t and x, t, respectively, identify a point of 
departure on the (lower) cone base line and the corresponding point 
of arrival at the cylinder surface. 

Since t^x is associated with t = 0, the first of the relations 
(36) gives 



cos 


max 


B 

B 


(37) 


which equation, as shown above, must be solved by a negative value 
(except with B = B where T max = 0) • The x- coordinate , xi = aa 1? of 
the point of contact along the line t = 0 can be obtained from the 
second equation (36): 


(x x - x^ ) sin a 


C sin t cos OC . 
max 


The right side here is not negative, so that x x 1 x^, as one expects. 

The cast shadow line on the lower portion of the cylinder thus 
begins at t = 0, a = and continues toward the nethermost generatrix, 
t = - Jt/2, which value, by the first relation (36), corresponds to 
t = - Jt/2. Thus, the abscisa, x 2 = a of the line's low point is 
given by 


(x 2 - x^) sin oc = (C - C) cos a. 

The right side again is not negative (since we assume that C ^ C) , and 
therefore x 2 = xb* 
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If either one of the differences (x x - x b ), (x 2 - x^) is larger 
than the cylinder length h = x c - x^, either t = 0 or V = - it/ 2 
cannot be reached by the cast shadow boundary on the cylinder, and 
the interval - it/2 ^ t ^ 0 must be curtailed accordingly. This can 
be decided with the aid of the system (36) which in fact is the equa- 
tion of the shadow curve s c in terms of the parameter T. One will 
have to determine at what value (or values) of t the line s c leaves 
the cylinder surface. If it does so, the curve s c (shown in Figure 3) 
would intersect with the line a = a c , thus extending the cast shadow 
region and removing a certain domain of negative T-values from the 
integration interval . 

For simplicity we assume that this is not the case, i.e., that the 
cylinder is long enough to have an underside area near its rear end 
fully "illuminated." Regarding the integration over cr, the upper limit 
will then be the natural limit a = a c . The lower limit, however, will 
vary according to the variation of cr along the curve s c whose equation, 
by eliminating r from the set (36), can be put into the form 


a = 


* +± 

b a 


C sin t + C 


•JT 


(B/B) 2 cos 2 t 


cotg ct = \|t(t) . 


(38) 


The values of a at the interval terminals, t = 0 and t = - it/ 2, are apt 
to be in any ratio depending on the axis ratios B/B and C/C. As one 
would also expect, there is an extremum at t = - it/ 2, either a maximum 
or a minimum. If (B/B) 2 ^ (C/C), a second extremum (a minimum) may 
exist. Figure 3 is a schematical sketch of the geometry involved, 
drawn up in the (x,a)-plane. Half of the natural r-interval must be 
excluded on account of the tangential shadow line, s t , while in the 
rest area the natural cr-interval is trimmed down by the cast shadow 
boundary, s c . 

The integration with regard to a proceeds from a - ^(?) to a = O c 
(increasing values). The second-order differential (34) then yields 
the first-order differential 


dZ = - 4q a B 3 sin 2 oi S — ~ T [cr - \|r(x)] dx. 

B 2 sin 2 x + C 2 cos 2 x C 

The quadrature in - n/2 < x < 0 will be carried out for two types 
of cylinders where it is relatively easy. 
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Cylinder Rear 


Illuminated Area 


Shadow Line s t ,T=0 


Shielded Top of 
Half 'Cylinder " 


= *(r) 


Cast Shadow Region-*)^ Cylinder Root 


Cone Base 
Cone Tip 


Figure 3. Shadow Geometry on Attached Cylinder (Schematic) 


A. The cylinders touch the major axis terminals of the cone base 
(B - B^. The curve s c descends here from the maximum to a minimum 

at t = 0. The functions and Z take the form 


+ — (C - C) sin r cotg a 


, _ . p_ h (* 1 A 1 _l+vl-A 

L 1 - A v 2 


f (C - C) 2 + — cot S «1 


( 39 ) 




where 


A = 



It is required that 

h ^ (C - C) cotg a; 


the cylinder must have sufficient length to accommodate the entire shadow 
boundary s c . Otherwise, Z assumes a different form. 

Two special cases stand out here: 

(a) C = C: the cylinder root is flush with the cone base. 

The curve s c reduces to o = cr^, meaning that there is no 
shadow cast on the cylinder. The expression (39) for Z 
loses the subtractive second term in the brackets (as it 
would in cross-flow, a = 90°). 

(b) C = 0: the cylinder degenerates into a horizontal plate 

of width 2B and length h ^ C cotg a. The cast shadow 
curve (36) is the ellipse 

* <>/»■ - 1 . 


A limiting process required as A -» 0 carries expression (39) into 

Z = Bq^ sin 2 a (4h - jtC cotg a) ^ BCq^ sin 2 OC 

For comparison, the Z- component for a fully impacted circular cone 
(fii ^ co, q = 0) follows from the second expression in the system (29) as 

Z = riBCq^ sin lot cos 2 w, (B = C), 

so that the lift of the composite body will be noticeably larger than that 
of the cone alone. The drag is less affected, although enlarged, too. 
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B. The cylinders are circular, C = B = R, BSCSR. The quantity 


is at most equal to unity, so that the square root in ^(t) is real. 
Expression (38) tells that 


Ox ~ t(°) = % + a J 1 " P 2 cot 8 a 

C “ R 

a 2 = \fr(-rc/ 2) = a, + cotg a . 

D Si 

The curve s c ascends from the minimum a 2 to the ordinate at t = 0 
(which is not a maximum). The rear end of the cylinder is fully illumin- 
ated if 


a 


c 


^ a 


15 


or if the cylinder length, h = a(a c - a^), satisfies the condition 


h ^ C n/ 1 - f3 2 cotg a, 


on which we will base the force computation. 

The differential dZ simplifies, since here 


B 2 sin 2 i + C 2 cos = R 2 ; 


one obtains from it 


Z = 4q R s in 2 QI -{-tt h + C cotg a 


16 C \8 2 


3?t R . fl fi 2 \ arc sin 3 \fl - 3 2 A- . 3 2 
1/1 77+ 0* o“ J ^ ' r ' 77 + J =T m 


J P 


V8 + 4 


4q R sin 2 a 
^00 


h + C cotg a 


16 C 3 + 15 + 160 + ' 


]} 
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The force decreases with decreasing cylinder radius and cylinder length 
(which, however, must not sink below the above lower bound; otherwise, 

2* would be a different expression). If the normal force coefficient is 
referred to the longitudinal cross section, 2Rh, 


C 


Z 


4 

3 


with (X = 90°, which is the known value for a circular cylinder in 
Newtonian cross flow (Ref. 2, p. 10). 

In calculating the moment, M, of the forces d 2 Z one best takes 
the cone tip (the origin) as the point of reference again. The integra- 
tion over o of expression (15) extends from a = ^(t) to a = a c * One 
finds that 


M = 2j q a 2 B 3 sin^a 




sin 3 r 


B sin^T 4- C cos 


dx. 


- jt/2 


The remaining integrand is elliptic except in the case A (B = B) where the 
square root in the expression for \Jr vanishes. Taking this case for 
illustration, 


p 

M = ~2i q ^ 2 ^ s in 2 a 

— J oo 





In 


1 -f 


+ 


° b a 


(2 - 3e 2 ) - 2(1 

— 


jfWf 


+ A 2 


3 - 5e 2 

3I 5 


a - g 2 ) 2 
21 s 


In 


1 + € 
i - e 


where 
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Since the minor cylinder semi-axis, C, is always smaller than the major 
semi-axis, B, the numerical eccentricity, e, of the cylinder cross- 
section can never become zero. The largest value that C may attain is 
P = C when the cylinder and cone bases are congruent (e = e) and M 
assumes the much shorter form 


51 


- M** «*;-■£) (i - i 
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In 


1-6 


which we will use in the following. The condition M = M* then gives 


x + x, 

C D 

x = o 


(40) 


Thus, the centroid of the cylinder is at its midpoint, as one should 
expect, since a cast shadow does no longer exist. 

The cylinder and cone moments both refer to the origin and can be 
vectorially added to obtain the total moment. The centroid of the total 
force can then be shown to have the abscissa 


x 


tot 


x*Z + x*Z 

9 

z + z 


when sought out on the body axis. If, again for simplicity, we assume 
the cone as fully illuminated (a £ go) , then from expressions (29), (31), 
(39) and (40) 


x 


* 


2it „ 1 - \l 1 - 

Cx, cos a 2 — 

3 D COSTCO 


X 




*b 2 


2 


sin OC 1 - 


1 - e g 
26 


In 


1+6 

1-6 


tot 


it C cos a, (1 - n/T - q) + (x^ - x^) sin a ( 1 - 


1 - 6 g 
26 


In 


1 + £ 
1-6 


where the abbreviations 


q = e 2 cos 2 co and A = (C/B) 2 = 1 - e 2 

have been used. The angle a can become 0°, but not 90°. In the first 
case x* ot = x* (cone centroid), since the cylinder, though flush with 
the cone, is not impinged at all. With very long cylinders, x* Qt £*. 
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If the cross sections are circular (e -» 0) , 


* 

*tot 


2 it Cx^ cos a + (x^ - x^) Sin a 


it C cos a cos 2 co + (x fi - x^) sin 


a 


Here, C is the common radius of the cylinder root and cone base. 


V. THE BIPARABOLIC CONOID 


With an elliptical cone the expressions (28) for the force com- 
ponents (a > oj) are lengthy. More concise formulations can be presum- 
ably achieved in the same flow if, in the (x,y) -plane, a sharp edge 
exists on the surface, preventing the formation of a tangential shadow 
line. For example, a body may be constructed whose cross sections 
parallel to the plane x = 0 (Figure 1) are bounded by two symmetric 
parabolic arcs facing each other and intersecting in the ground plane 
z = 0. Consider the two parabolas (parameter a > 0) 


y 2 = 2a (z Q + z) 


in the plane x = x D > 0, with the upper (lower) sign applying in the 
upper (lower) half-plane. The two finite arcs emerging with the restric- 
tion OS |z| S z Q enclose the area 


8 

3 


\I2slZo, 


since z = 0, y = ± -J 2a z Q are their common points, and z Q is their 
height. Writing now 

y 2 = 2a (x tan co + z) (4-1) 


we obtain a body whose cross-sectional areas are biparabolic and taper 
off toward zero if x -» 0; there is a tip at the origin. The body roughly 
resembles the elliptical cone illustrated in Figure 1, with which it will 
be compared. The plane y = 0 intersects with it in the straight lines 

x tan w + z = 0 
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which have the slope tan w and, incidentally, are the only straight 
lines on the surface. The planform (in the plane z = 0) is the rim 
parabola 

y 2 = 2ax tan u; (42) 

it is not a triangle as with the cone. The tip angle, go, the body length, 
Xfo, and the parabola parameter determine the length of the ,f semi-axes” 
(Figure 4) 



Figure 4. The Biparabolic Conoid 
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If B and C remain small while xb grows large, the body assumes an arrow- 
like appearance. If, on the other hand, the parabola parameter (a) is 
very large and co remains small, a wing-like structure emerges. The volume 
of the body is 


16 

15 


B C 


v 


it 

almost equal to that of an elliptical cone (j BC x^) of equal length and 

equal base axes. The bigger planform is offset by the smaller cross 
sections (their area ratiosare 4/3 and 8/3jt, respectively). 

One verifies readily that the surface (41) can be given the point- 
wise representation 



(43) 


The natural range of the variable a is 


0 ^ a ^ a b . 


(44a) 


corresponding to 0 ^ x ^ x^. Again, from 


0 ^ ± z ^ x tan co, 
tan co = = 0 , 


and thus. 


- \l tan co ^ t = \l tan co . (44b) 


The negative values of t cause y to be negative, while the parenthesis in 
the expression for z is always positive (or zero). 
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Both sets of the parametric curves (10) are planar, since the 
coordinates (43) depend in the second degree only on a and t. The curves 
o - o * (second family) are the biparabolic arcs. The curves t = t* con- 
tain the origin (with a = 0) and are therefore in the planes 


z 


± (tan go - 



y = 2acnr* (arbitrary) 


through the y-axis. The angle, X, they make with the (x,y)-plane is 
given by 


tan X = ± (tan go 


■), 


which relation, since t* may be any one value within the range (44b) imparts 
a geometric meaning to the variable r. 

For finding the shape of the curves t* = const., consider one of the 
planes, fixed by, say, a positive value of X (upper sign). The rotation 


x = £ cos X - £ sin X 

y = Tj 

z = £sinX + £ cos X 


carries this plane into the plane £ = 0. By the system (43), the curve 
t* in it is described by 


£ cos X = 2aa 2 

T) = 2aax*; 

elimination of a gives 

T| 2 = 2a t* 2 £ cos X = 2a (tan go - tan X) £ cos X. 
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The first family therefore consists of parabolas cut out of the conoid by 
planes through the y-axis . The planform (42) is one of them, correspond- 
ing to X= 0, or to x* = tan go. It will be observed that, quite generally, 
the right parabola branch (y > 0) is associated with a positive value of 
x*, the left one with the oppositely equal negative value. Moreover, the 
same two values appear with two parabolas, one on the upper half (X > 0) 
of the body, the other on the lower half (X < 0) . The force components 
must be set up separately for the two halves (the body is composite, in 
fact) . 


The straight ridge lines are characterized by X = ± go, or x* = 0. 
The parabolas grow more slender if X approaches ±go and finally degenerate 
into a straight line. At the same time | x*| becomes smaller and is zero 
at the end. 

From relations (12) the interior normals are found as 


rin = 


sin go 


0 1 + (2or cos go)' 


2ax cos go 


Jl + (2crx cos go)' 


n_ = + 


cos to 


\l 1 + (2ax cos go ) 2 


while 


q 2 2 

N - ~ a n/I. + (2ax cos go) 2 . 
cos GO 

The upper sign refers to the body part in the upper half space. 

If the flow vector, as with the cone, is given by expression (25) 
(the particles arriving from below and behind parallel to the (z,x)~ 
plane), equation (18) assumes the form 


sin (go + OC ) = 0. 

There can be no tangential shadow line on the lower part of the body (as 
is evident from inspection). On the upper part, a shadow line s c exists 
with a = go only; it coincides with the upper ridge-line, encloses no 
area and does not affect the forces (C p = 0 along it). 
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The shadow cylinder through the rim parabola (42) is found to have 
the rim points only in common with the body so that a cast shadow curve, 
in the proper sense, does not exist either. However, the rim evidently 
will act as shadow boundary in certain circumstances. To discuss these, 
we consider a particle path piercing the ground plane z = 0 on or outside 
the rim parabola. It intersects with the upper surface, i.e., the moving 
particle will strike it, if 


tan oo cotg a, ^ 1. 


The proof, omitted here, starts out with a point on the upper surface and 
follows the path backward toward its intersection with the ground plane. 
The choice of the point is completely arbitrary so that we can say: as 

long as a i oj, the entire upper surface will be impinged upon; with OC > oo 
none of it will. Instead of the a-dependent shadow line T t (cone formula 
(27)), we have here the fixed boundary 


x 


± \l tan oo 


with OC > oo. The force expressions will be simpler accordingly. Applica- 
tion of formula (14) gives 


lbq^ a 2 sin 2 (oo + a) 

d 2 P = - — ; — 7 — (i sin oo - 2iaT cos oo + k cos oo)dadT 

— cos oo 1 + 4a r cos oo J 

where the upper sign refers to the elemental forces attacking the body's 
upper half. On integrating over the t- interval 

-s/ tan oo ^ r < \J tan oo 


the j-component is found to vanish (dY = 0) leaving us with 

16q a 2 sin 2 (oo + O') , 

dP = 5 cj arc tan (csv2 sin oo) (i sin oo + k cos oo) da. 

— cos^oo 
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Further integration over the natural range 0 ^ a ^ results in 


X = 2q a 2 

^■oo 


sin 


2 (C0 + a) 

— 5 

COSTCO 


[(1 + m 2 ) arc tan m - m] 


Z = + X cotg co 


where 


m = a, n/ 2 sin 2co = / — sin 2co . 
b v a 

With the use of this abbreviation, the conoid's base area, to which the 
force component will be referred, may be written as 

— a 2 (m/cos co) 3 . 

Take at first a > co (lower signs apply alone). Then, by the system (30), 


D 


3 sin 3 (co + a) 

° sin co 


= 3 sin 2 (co + a) cos (co + a) 
C L 2 A sin co 


a > co 


where 




+ m 2 ) 


m 


arc 

3 


tan m 


m 


These expressions are considerably more compact than the corresponding 
cone formulas arising with the systems (28) and (30). It is seen that 
with a sufficiently large angle of attack (co + OC > 90°) the lift coef- 
ficient is negative indicating that the aerodynamic force then seeks to 
drive the body beneath the instantaneous flight direction. 
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With a ^ to, the upper part of the body is impinged upon also, and 
the force components involving the upper signs must be added to the lower 
part components: 


C * 3A cos a [sin 2 to cos 2 a + 3 cos 2 co sin 2 a] 

D 

CL = 3A sin a [cos 2 to (cos 2 a - sin 2 a) + cos 2 a (cos 2 co 
JL 

The lift coefficient here is positive up to at least to = 45°, with ot -» 0° 
up to to = 54°. The corresponding figures for the elliptical cones are 
to = 36 0 23' and to = 4l°49 f , respectively. 

With ot = to both sets of expressions give 



= 12A sin 2 to cos 3 to 
C = 3A cos to sin 2to cos 

j_i 



a = to. 


The force coefficients thus remain continuous if the upper surface passes 
from full illumination (fit < to) into full darkness (O' > to) . Indeed, with 
a = to, X U pper = z upper = meaning that the flow vector, is then 
tangential to the upper surface everywhere. 

One may compare the force coefficients to those of an equivalent 
cone having the same base area, the same length, and the same tip angle, to. 
In these circumstances the cone base major semi-axis, B, is smaller than 

B (Figure 4); B = B, while C = C. As a consequence, the volume of the 

conoid is by 20 percent larger than that of the equivalent cone. 

Using the formulas (29) one finds that the cone's lift decreases 
with its length, x^; it remains non-negative even with x^ 0, if 


tan 2 a ^ 


2 - 3 sin to 
1 - sin to 


(This must be satisfied in addition to tan Ot i tan to.) 


& 
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If and Rq denote the lift and drag ratios, respectively, of the 
conoid and cone, one can, for small angles of attack, derive the expres- 
sions 


— — AK (2 cos 2 to - sin 2 u))(l + K sin co ) 
L ^ 2 K(2 cos to - sin 2 to) - sin w 

3 

R_ -- AK sin to 
D 2 

where 

K = n/I + (x b /B) a . 


Use has been made here of the relations 


q = cos 2 o) - (3jt/16) 2 m 2 = 1 - K 2 sin 2 to 


which follow from the equivalence as stipulated above. 

As long as the ratio R^ is positive, it increases with to, ranging 
from Rl ~ 1 (when to = 0) to R^ » oo (when the cone lift is zero). 

The drag ratio increases from Rq ~ 0 (at (0 = 0) to R D « 1 (at (o = jt/2). 
Thus, if a « 0, the conoid’s lift is larger than that of the equivalent 
cone; its drag is smaller, although the conoid has the larger volume. 

Comparisons when OC is appreciable, or when a > to, have not been made; 
the formulations grow rather unwieldy then. 

The location of the centroid on the conoid's axis again does not vary 
with angle of attack: 


X* = 


*b 


r m + j m 3 - (m 2 +1) arc tgm' 

<1 + (1 - 2 sin 2 to) —27 ; — 

1 v 7 m [-m -f- (m + 1) arc tgm] 


2 cos^to 
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However, there is a dependency on the base shape not found with the 
elliptical cone. Note that the parameter m may be written as 

2x b . 

m = a— s m go. 


while the base area is 


8 a 

3 Bx b 


tan go. 


If, on varying go at constant body length, this area is kept constant 
rather than the semi-axis B (B = B cotg go, B = const., C = x^ tan go), 
the body volume 


— K v 2 

15 B ^5 


tan go 


will also remain constant, but the general appearance will change from 
wide and flat (wing-like) at go « 0 to narrow and high (rafter-like) at 
large values of go. It can be shown that in these circumstances the 
centroid abscissa keeps increasing with increasing values of a), just 
as it does with the elliptical cone. 

With two special tip angles the location of the centroid is inde- 
pendent of the choice of B. 

A limiting process shows that with go = 0 (m = 0) 



which is somewhat smaller than the cone value x ,f = j x^. The center of 
gravity is at x = y x^. Like the elliptical cone, the conoid is not stable 


if it is wing- like. One can calculate that stability requires go > 28°, 
while with the cones the tip angle must exceed ~ 20°. If go = 45°, x* = x^. 
The centroid then coincides with the center of the base. For the elliptical 
cones this occurs at go ^ 35°. It appears that the outward movement of the 
centroid with increasing tip angle is slower with the conoid than it is 
with the cone. 
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VI. A BLUNT BODY 


The lemniscatic body reported on in Reference 3 was investigated 
for three reasons; first, to test the general method with a surface of 
more intricate description, second, to show that the tangential shadow 
boundary, if not coinciding with a surface coordinate line, may call for 
close attention as to integration limits and that some of the quadratures 
required may have to be done numerically, and third, to study a body some- 
what similar to the Apollo capsule. 

The latter purpose is served better when we now consider a body com- 
posed of a right circular cone and a spherical base cap whose radius, R, 
is rather considerably larger than that of the cone base (B) . The mathe- 
matical formulations to follow are valid down to R = B, the cap becoming 
a half-sphere. The angles p and p identified in Figure 5 determine the 
coordinates of the point P on the sphere: 


x = R (1 - cos jj,) 
y = R sin [i cos p 
z = R sin [x sin p. 



38 


Figure 5. The Blunt Body 


They stand for the variables o and T in general theory. The parametric 
curves p = p* and \i = \x* are the meridians and the circles of constant 
latitude on the sphere, respectively. One sees easily that the condition 
(11), A = 0, for planar curves is satisfied. Let \i = p. b denote the lati- 
tude of the cone base. The natural intervals then are 


0 ^ H ^ ^ 

0 ^ p ^ 2jt. 


Expressions (12) give 


n x = cos \x 9 n 2 = - sin pcos p, n 3 = - sin |i sin p 
N = R 2 sin |a. 


If the particles arrive from the lower left. 


v_ = i cos a, + j sin a. 


and the condition (18) for the existence of a tangential shadow line on 
the sphere assumes the form 

cos p = cotg a cotg p, t# 

One sees from this that the cap will be fully illuminated if 


Otherwise, the force and moment integrands will contain a term 
arc cos (cotg OC cotg \x) , equally unpleasant as the term arc cos (cotg OL 
cotg 3a) that can never be avoided with the lemniscatic body (Reference 3). 
It is likely to force a numerical procedure upon us when integrating over 
jj.. Fortunately, the Apollo cap is very flat (R » B) so that will be 
a small angle. Since large angles of attack can be excluded, the tangen- 
tial shadow line on the cap may be taken as non-existent. The force dif- 
ferential (14) can then be integrated without difficulty, since the fixed 
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natural limits apply. If the force coefficients are referred to the cone 
base area, one finds that 

C sin 2 a sin 2 p^ + (1 + cos 2 ^) cos 2 OJ 

C Y = sin a cos a sin 2 p^ 



Figure 5 indicates that the unit "a" treated as arbitrary in Sec- 
tions III and IV is taken here as the cone length. Since x = p for a 
circular cone (Figure 2), its point-wise representation may be given as 


x = a (o + 1 - a) 
b 

y = B a cos p 
z = B a sin p. 


This system follows from the former system (21) by appropriate coordinate 
transformations . For greater clarity, p has been written for p, although 
the meaning of the angle is virtually the same as that of the angle p used 
with the sphere. The variables a and p range in 


0 ^ a ^ 1 
0 ^ p ^ 2 xc . 


The values a = 0 and a = 1 identify the cone tip and the cone base, 
respectively, so that, when integrating with respect to o, one moves 
from tip to base. 

From expressions (12), with the use of B = a tan gj, 
n x = - sin (a), n 2 = - cos p cos co, n 3 = - sin p cos go 


s in go 
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Condition (18) for tangential incidence becomes 


cos (5 


tan co 
tan a 


An illuminated area on the cone can exist only if the angle of attack is 
at least equal to the tip angle (O' ^ co), as one would also infer from 
inspection. If p t = p* is one solution of this equation, p t = 2^1 - p* 
will be the other one. 

The shadow cylinder through the base periphery can be shown to inter- 
sect with the cone in a line wholly confined to the tangential shadow 
area, so that the integration interval is not curtailed any farther by a 
cast shadow boundary. 

The force differential (14) assumes the form 


d 2 P = - 2q 


B 2 a 


(sin oo cos a + sin PC cos oo cos p) £ 
sin oj 


(i sin oo 4* j cos co cos p + 


+ k cos co sin p) da dp. 


It needs to be integrated only if a ^ oo (otherwise, the cone does not con- 
tribute to the force). The ranges 0 ^ a ^ 1 and p* ^ p ^ 2 it - p* must be 
used. The coefficients, referred to the area rdBf, may be written as 


c x - 


sin 2 q cos 2 gj 


(1 + 2 cos 2 p*) (it - p*) + 3 sin p* cos p* 


C Y = 


2 sin 2 q cos 2 co 




(tc - p*) cos p* + S3 "^ - P (2 + c‘os 2 p*) 


cotg oo 


C 2 * °- 


Tt 

The brackets here are found to be non-negative in the interval, < 2i 
allotted to p*; the first varies from at/2 to zero, the second from 2/3 to 
zero. Thus, if the cone is struck by the particle stream, it will sub- 
tract from the chordwise forces and add to the normal forces of the cap 
alone. 
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The centroid of the cap is at x* = R (even if a tangential shadow 
line should exist on it), that of the cone is at 

** " *b + (‘ - 1 - 

as follows from expression (31) allowing for the changed orientation and 
denotation (former x^ = present a). Since, from Figure 5, 


X^ = R(1 - COS (jL^) 


the arm of the total moment will be 


x 


tot 


RY + x*Y 
Y + Y 


a 

R + - 


R cos a. 

b 


2 a 

3 cos 2 co 


Y + Y 


Y 


= R, if there are no forces on the cone. 


Since 


I < ± 

Y = 3Jt 


tan a cos 2 go sin 2 ^, 


the centroid of the Apollo configuration, where jjl^ is a small angle and 
a < R,will never stray far from the cap f s center. On the other hand, 
the center of gravity will be near the cone base. Since x^ « R, the 
configuration should exhibit strong stability in the Newtonian flow 
regime. 

The same can be expected of the Gemini capsule where the rear end of 
the cone is replaced by a coaxial cylindrical piece. To be sure, if the 
conical middle part is struck by the particles (a > w) , the deflected 
stream will come in contact with the cylinder surface, and the Newtonian 
analysis cannot be applied with confidence. However, the presence of the 
cylinder should tend to drive the center of pressure still more backward 
and thus to further stability. 
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VII. THE MODIFIED APPROACH 


With some bodies of plane or axial symmetry and with the circular 
cylinder in symmetric cross-flow, the Newtonian results have been shown 
to improve if one sees to it that the pressure coefficient assumes the 
exact value at the stagnation point where it is usually (because rela- 
tively easily computed for O' = 0) . It may be expected that the expres- 
sion thus gained will also hold good for small angles of attack. At 
least one corroboration of this surmise exists in the pressure distribu- 
tion around a circular cone (go = 10°, QC = 6.7°) where the modification 
amounts to a 4 percent increase in values that are already satisfactory 
on the whole when computed from shockless impact theory. 

With the overall angle of attack zero, the angle, go, at the stagna- 
nt 

tion point will be the local angle of attack so that w 2 “ a stag* ^ 
modified formula then will be written as 


C = c * = c * 

p p sin^co p 


(v • n ) 2 
sin 2 to 5 


(45) 


where Cp is the pressure coefficient at the stagnation point. If 
a} = o£ ta g, Cp = C* as desired. The value of Cp can be calculated from 
shock transition relations and depends on the ratio of specific heats 
( 7 ) in the gas and on the Mach number, M^, of the undisturbed flow. 

In the case of blunt bodies (for which expression (45) was first 
suggested in Reference 4) sin go = 1 and 


C* 

P 


2 1 \ (y + 1 

7 Lv 2 




+ 1 V -1 

- 7 + v 



(46) 


With infinite Mach number in a diatomic gas, C* = 1.84, which figure then 
replaces the factor 2 of simple impact theory. The values decrease with 
decreasing Mach number (C* =1.64 for M^ = 2) , at first very slowly; in 
the hypersonic region M^ > 6 the figure 1.84 may be used throughout with 
a small error in the second decimal place (7 = 1.4). Very satisfactory 
results have been obtained regarding the sphere, ellipsoid- arid sphere- 
capped circular cylinders, and a sphere-blunted circular cone; they were 
somewhat less accurate with the cylinder in cross-flow. In all cases, 
however, they surpassed those obtained by another method (Busemann*s pres- 
sure relief approach). 
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With plane symmetric bodies having a sharp leading edge to which 
the shock is attached, one may use the zero incidence stagnation pressure 
of the wedge which, although it cannot in general be written down explic- 
itly, assumes a convenient form when the cosine of the shock-body angle 
is sufficiently close to unity. Expression (45) then emerges as 


C 


P 




+ 


4 

M^ sin^co 

00 


(v • n) 2 . 


(47) 


If the Mach number approaches infinity, the bracket approaches (y + 1); 
the factor 2 is then replaced by 2.4 in a diatomic gas. 

With y = 1.4, the formula (47) worked well and better than Busemann 1 s 
method for the wedge itself and for a symmetrical pointed airfoil profile. 
With the latter and y = 1.05, however, the modified Newtonian formula gave 
consistently too high pressure values and was inferior to the pressure 
relief approach (which resulted in figures somewhat too small). 

The surface of a pointed body of revolution may, near the tip, be 

approximated by that of a circular cone with the same half opening angle 

co. The latter 1 s relation to the angle, a g , of the attached shock is 

involved. As a rule, numerical calculations are necessary, unless both 

co and a are small. In this case the approximate expression 
s 


c* 


+ 1 


- 


is derived in Reference 1 
K s = Moo c s being given as 


1) + 2(K - K) 2 (48) 

(7 - 1) + -gz 

s 

(p. 116-118)'", the relationship of K = and 


K . t 

_s = 7 + 1 + 

K 7 + 3 


(49) 


Much of the factual information assembled in the present section is also 
taken from Chernyi 1 s book. 
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2.08 


If M co 


-» °°» 


C p (y ± i) (y ± 7) 
w 2 (7 + 3) 2 


with 7 = 1.4. The excess over 2 is markedly less than in two dimensions. 
For the circular cone itself and 7 = 1.405, the approximation of Cp is 
very good up to go = 20°, 30°, and 40° if K 1 2, ^3, = «>. it breaks 
rapidly down for K < 2, the error amounting to -8 percent at K = 1 and 
0 ) = 5°. Expression (48) offers an equally satisfactory approximation of 
the ratio o s /u> in terms of K; with co up to 10°, it is close even with 


K = 1. 


A check was also made with an axisymmetr ic ogive (gj = 16.26°, Moo “ 8 , 
7 = 1.4). The zero incidence meridional pressure distributions as com- 
puted from Newton's modified formula and from the (more exact) numerical 
method of characteristics were practically identical. 

For bodies like the elliptical cone and the biparabolic conoid which 
are not of rotational symmetry,, the modification of the factor 2 must be 
judged on the basis of the wedge and circular cone results. The flatter 
these more irregular bodies become at a given value of w, the more one 
may be inclined to cautiously upgrade the relative low cone correction. 

The blunt Apollo-like configuration induces no uncertainty; the modified 
pressure coefficient will here be smaller in accordance with the general 
expression (46). 
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